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Photon Absorption on a Neutron-Proton pair in 3 He and 4 He 
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We have developed a phenomenological quasi-deuteron model to describe photodisinte- 
gration reactions of 3 He and 4 He at intermediate photon energies, and compared the results 
with experimental data obtained by TAGX group in which contributions from photon ab- 
sorption by pn pair were selectively observed. The data can be reproduced reasonably well, 
however there still remain certain discrepancies between calculated results and data. 
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\Q ■ §1. Introduction 
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Photodisintegration reactions of light nuclei have been studied for a long time 
both theoretically and experimentally. Experimental results of these reactions are 
naively expected to be reproduced precisely by theoretical works since the reactions 
are induced by the well-known electromagnetic interaction and the structure of few 
nucleon system is believed to be calculable. Despite of this expectation, there still 
remain certain discrepancies between experimental data and theoretical results. We 
are interested in the origins of these discrepancies. 

The photodisintegration of nucleus was studied theoretically using the quasi- 
deuteron model by Levinger& and later by Futami and MiyazimaEf. They expressed 
the photodisintegration cross section of nucleus using those of deuteron with a certain 
factor, so called Levinger's factor, which accounts essentially for two effects, the 
relative wavefunction of pn pair and the effective number of deuteron in the target 
nucleus. Recently, the Levinger's factor was updated and determined for many nuclei 
using the modern data of rms radius a> , which enables us to obtain information of the 
relative wavefunction of pn pair. The accuracy of the model was also investigated for 
wide range of momentum transfer for trinucleon systemic . These results reproduced 
the experimental data qualitatively well indicating that the photon absorption by 
pn pair is important in photonuclear reactions. 

On the other hand, the deuteron photodisintegration reaction, which is treated 
as the elementary process in the quasi-deuteron model mentioned above, was studied 
in microscopic way by Leidermann and ArenhovelEf. They calculated the two-body 
photo break-up of the deuteron above pion threshold including A degrees of freedom 
explicitly. Homma and Tezuka also developed another microscopic model of the 
deuteron photodisintegration reactionB 1 . These models described experimental data 
well at intermediate energy region, E 7 > 150 ~200MeV. 

The microscopic models were applied to J/y, p) reaction for several nuclear tar- 
gets assuming the quasi-deuteron picture B*' IIP. In the applications, they assumed 
two nucleon photon absorption processes expressed microscopically, and obtained 
the spectra of ejected protons after absorbing the photon, including the distortion 
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effects as the attenuation factor. They did not include the contributions from neither 
spectator protons nor ejected protons after multiple scatterings which are expected 
to contribute in lower momentum region of proton spectra. Thus, their models could 
only be applied to emitted protons with higher momenta. 

Since 1987, new experimental data of the jphotodisintegration reactions of the 
light nuclei have been obtained by TAGX group W . They have obtained kinematically 
complete data of 3 He(7,pn)p and 4 He(7,pn)d for the first time&'Br. They have 
striven to separate the whole events mto two- and multi- nucleon absorption processes 
by the momentum ordering method c3. We are very interested in the data since the 
contributions only from the photon absorption by pn pair are selectively observed 
which should be calculated precisely in the quasi-deuteron picture, and since the 
data are completely exclusive which do not include any extra particles like pions in 
the final state. Thus, we think it extremely important to investigate this data to 
know whether it can be understood theoretically or not. 

So far, the 3 He(7,pn)p reaction was studied theoretically by Wilhelm et aZ.ll^P'll^P . 
They considered both triplet (i.e. deuteron channel) and singlet configurations for 
initial pn pair in the nucleus, and included photon interaction with nucleonic, mesic, 
and A currents. They calculated the cross sections and spin observables, and found 
that the data of total cross sections deviated from the theoretical results at lower 
energies (E 7 < 200MeV). Since their results for differential cross sections were not 
in the form which could be compared with the data directly, there have been no 
comparisons between data mrd theory for angular dependent observables yet Ed). 

The 4 He(7,pn)d datatH' show that the total cross section has strong energy 
dependence around E 7 = 150MeV which have not been investigated theoretically. 
Theoretical results by TezukaLP are consistent to the data for higher energies but 
there are no theoretical calculation which can be compared to the data in whole 
incident photon energy region. Their results include the attenuation factor to account 
for the distortion effects and the final nuclear state is not specified to deuteron. 

In this paper we attempt to understand these reactions, 3 He(7,pn)p spectator and 
4 He(7,pn)d<j PeC i a to r , theoretically. For this purpose we investigate both 3 He and 
4 He photodisintegraion reactions within a same theoretical framework and try to 
reproduce all data of the both reactions simultaneously. We calculate all observables 
in suitable form which can be compared to the data directly. Since whole events 
are classfied carefully by kinematical information of all final particles in the TAGX 
experiments, we could also have possibilities to check their classifcation scheme of 
the data independently by studying the theoretical results assuming the two body 
photoabsorption. The present work using the phenomenological model is expected 
to give important guidelines to fully microscopic studies which should be done as a 
next step. 

We describe our model in the section 2. Numerical results are compared with 
the data in section 3. We summarize this paper in section 4. 
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§2. Formalism 



In this section we describe our phenomenological model which is applied to both 
3 He(7,pn)p<j perfator and 4 He(7,pn)d sperfator reactions. Since we would like to calculate 
all observables in suitable forms to compare with the data directly, we evaluate the 
final-state three-body phase space exactly. The cross section for A(7,pn)B reaction 
can be written as; 



da = X1/2 ^ Q)M%) I Ta I 2 5\ PB +p P + p n - PA - Pl ) 

M B M p M n M A 3 3 3 



EbE p E„ 



(FpBcPppdrpn, (2-1) 



where A indicates the target nucleus, 3 He and/or 4 He, and B the spectator particle 
in the final-state, p sp ectator and/or d spec t a tor- The s corresponds to the Mandelstam 
variable in the initial 7+A system. We integrate the final-state phase space properly 
for each experimental result. 

The A(- • •) is the Kallen function, which is introduced to normalize the cross 
section properly by the initial photon flux, defined as; 

A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 26c - 2ac. (2-2) 

In our model, we consider only deuteron channel for the initial pn pair in the 
nucleus and express the square of the amplitude of photon absorption by the nucleon 
pair using the differential cross section (jj^j of d(7,p)n reaction as; 

I T hpn P n) |.= {%) 0, Ml^f^.^- y ( 2 .3) 

where all kinematical variables are evaluated in the center of mass frame. We use 
the phenomenological fit of the {jfj^J datatij), which reproduce the d(7, p)n cross 
section well in the wide energy region, E 7 =20-440 MeV. 

In order to calculate the observables for He target cases, we postulate the quasi- 
deuteron picture and describe the square of the 7 absorption amplitude by He nucleus 
as; 

I T A | 2 = | T(jpn -> pn) \ 2 \ ^(j>B)pn-B | 2 , 

where we evaluate the | T(p/pn pn) \ 2 using the Mandelstam's s and t variables for 
two nucleon system which absorb the photon. The rp is the relative wave function 
between pn pair and spectator particle B in He nucleus, which is assumed to be the 
Gaussian form as; 

Hpb) P u-b = (2tt) 3 / 2 ' exp (-^f) , (2-4) 
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where pf means the Fermi momentum between pn pair and the spectator particle 
B and is roughly estimated to be 146 MeV for 3 He and 205 MeV for 4 He from the 
data of the charge radii. p_B is evaluated in the laboratory frame. 

Here, we would like to describe some features of the present model. The advan- 
tages of our model are that we can calculate any cross sections which can be compared 
to the data directly for both 3 He and He target cases and that we can calculate 
observables in whole energy region of interest using the experimental information of 
the deuteron photodisintegration reactions. However, we do not include neiter the 
effect of the deuteron spectroscopic factors nor the effects of the compact wave func- 
tion of the pn pair in the target He. Both effects are expected to be insensitive to the 
incident photon energy and to change the overall normalization of the cross sections. 
Thus, we feel free to introduce a constant factor to each reaction as a parameter to 
normalize the absolute value of the cross sections when it is necessary. Actually this 
factor is considered only in the case of 4 He(7,pn)d spec t a t OT . reaction as we will see 
later. For 3 He target cases, it is not necessary to consider this normalization factor 
and this could be accidental. In the model, we do not include the contributions from 
the singlet pn pair in the initial state, which were evaluated in refs.@'0) and shown 
to be around 50 fib at E 7 = 100 MeV for 3 He(7,pn)p sperfator total cross section. This 
contribution decreases monotonically as a function of E 7 and is known to give minor 
contribusion to the photodisintegration reactions at higher energies. We should also 
mention here that the final state interaction is partly included in the model from the 
beginning since we have used the fit to the data of deuteron case. The final state 
interaction effects for spectator particles are roughly considered in section 3.3. We 
compare the results with existing data and try to understand featuers of the data. 



§3. Numerical Results 



3.1. 3 He(j,pn)p spectat or reaction 

In this section we compare our calculated results with experimental data of the 
3 He(7,pn)p< jpectator reaction taken by TAGX groupt^'E-lP. 

First we consider double differential cross sections (DDCS) for emitted protons, 
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We expect to be able to fix the pf value defined in the last section by 
the widths of the DDCS. In Fig.l, theoretical results of DDCS are comapred with 
the data at three different angles for each incident photon energy. In all cases in 
Fig.l, experimental data show two distinct peaks which correspond to spectra due 
to the spectator protons and energetic protons paticipated in photon absorption 
process. The peaks at higher momenta have stronger angular dependence than the 
spectator contributions and become small rapidly at large angles. Our calculated 
results with pj=150MeV reproduce all DDCS data reasonably well without any 
adjustable parameters. In order to see the sensitivity of the calculated sepctra to 
the pf value, we also show the results with pj=200 MeV in Fig. 1. We see that the 
spectra are affected by the pf value as we expected. However, it seems difficult to 
determine the pf value precisely from these data. Thus, we use the p/=150MeV, 
which is consistent to the charge radius of 3 He, in the following calculations. 
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The differental cross sections for protons are shown in Fig. 2 for four incident 
photon energies. Except for the .E 7 =325 MeV case, the data are reasonably well 
reproduced again without any parameters. At £" 7 =325 MeV, the observed differen- 
tial cross section is small at forward angles in contrast to calculated results. The 
differential cross sections are insensitive to the pf value and are qualitatively the 
same for all pj values which are consistent to the DDCS data shown in Fig. 1. 

We show the total cross sections of the 3 He photodisintegration by two body 
(pn pair) photon absorption process in Fig. 3 as a function of the incident photon 
energy. Experimental total cross sections are obtained by integrating the differential 
cross section data in Fig. 2 over the observed kinematical region. The contribution 
from unobserved kinematical region is corrected for by extrapolation using the two- 
nucleon absorption modelUHP . Theoretical results only include the deuteron channel 
of initial pn pair in the target. As we can see from the figure, the data show the peak 
at E 7 =225MeV, which seems due to the A excitation, while the calculated results at 
275MeV. This difference could be explained by inclusion of the siglet configuration 
of pn pair in the initial state which was evaluated by Wilhelm et and 
shown to be larger for lower photon energies. However, the energy dependence of 
the theoretical results will become more different from the data by including the 
singlet contribution at .E 7 =125 ~ 205 MeV. In this energy region, the data increase 
monotonically with photon energy while the theoretical results may be flat or even 
have opposite energy dependence by including the singlet contribution. Our results 
are consistent to those calculated by Wilhelm et aZli^'EiP for initial pn pair in the 
deuteron channel. 

3.2. i He('j,pn)d spec tator reaction 

Total cross sections of the HeIs/,pn)d S p ec t a tor reaction were measured at several 
photon energies by TAGX group and shown to have strong energy dependence 
around £' 7 =150MeV. This energy dependence is almost like a 'discontinuity' to lower 



energy data taken by different groups 



y . We first calculate the total corss 



section with pj=200MeV, which is consistent to the charge radius of 4 He. We show 
the results in Fig. 4 and Fig. 5 together with the data. As can be seen in both 
figures, the calculated results are larger than the data at lower energies and smaller 
than the data at higher energies, and do not reproduce the strong energy dependence 
at 150MeV. 

Then, we vary the pf in wide range to investigate the possibility to reproduce the 
energy dependence by changing the wavefunction of 4 He. We also multiplied a factor 
to the calculated cross sections to reproduce the A peak. The results are shown in 
Fig. 6 and Fig. 7. We find that the our results do not reproduce the experimental 
energy dependence. It is difficult to reproduce the observed 'discontinuity'. 

We would like to mention here the recent work by Efros et a/.EP , in which the 
total 4 He photodisintegration cross section is studied in the giant resonance energy 
region. They have shown both the total and the two-body decay contributions, which 
are He(7,n) 3 He and 4 He(7,p) 3 H reactions, separately for 4 He photodisintegration 
cross section. We can estimate the absolute value of cross section of our case (3-body 
decay) by subtracting the 2-body decay contribution from the total. It is about 1 
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mb at E 7 ~ 40 MeV which is even larger than our largest result(~ 500 fib). If this 
is correct, the older data at lower energies may include large errors. 

3.3. Final state interaction effect 

In our quasi-deuteron model, since we have used the d(~/,p)n data to evaluate 
the amplitude of the photon absorption process by pn pair, the effects of the final 
state interaction between proton and neutron, which absorb the photon, is included. 
However, we have not included the distortion effects between the spectator and the 
other particles. We think it necessary to evaluate the final state interaction effects 
for spectator particles carefully for more realistic calculations. Since the two body 
photon absorption events are selected by kinematical information by TAGX group, it 
could be difficult to separate multi-nucleon photo-absorption events from two-body 
absorption events affected by the final state interaction, especially at low E 1 cases. If 
the spectator has relatively larger momentum because of the final state interaction, 
it could not be seen as the spectator. 

In this paper, we have estimated the final state interaction effects in a simple 
way to see sensitivities of cross sections by changing the ip pn -B defined in eq. (2.4) 
as follows; 



^'(PB)pn-B = / dr B [l-fexp(-r 2 B /rl)]e^(-ir B -p B )^(r B )pn-B, (3-1) 



where ip(r B ) pn _ B is the relative wavefunction between pn pair and spectator particle 
B in coordinate space, of which momentum space expression is given by eq. (2.4). 
Parameters / and ro are introduced to distort the plave wave of spectator particle, 
exp(— ir B ■ Pb), in the final state. Results without distortion effects correspond to 



We calcluate total cross sections for both 3 He and 4 He cases with changing the 
parameters in the range of — 1 < / < 1 and 0.5 < ro < 2.0 [fm]. We find that 
the absolute value of total cross sections changes significantly by this modification. 
However, in the present simple model, energy dependence of the cross sections is 
almost the same as those with / = case. As a natural expectation, the distorion 
effects have energy dependence. Thus, we need to include realistic energy dependence 
in the distortion effect for further studies. 



In this paper, we have investigated the photodisintegration reactions 3 He(7,pn)p s . 
and 4 He(7,pn)d sp which have been observed by TAGX group. We have investigated 
both reactions within a same theoretical model and comapre the calculated results 
with experimental data. In the model, we have assumed the quasi-deuteron mecha- 
nism to describe the photon absorption amplitude using the experimental informa- 
tion of deuteron photodisintegration. We have introduced the relative wave function 
between pn pair and spectator particle in the target He and treated the phase space 
integration carefully in order to calculate the observables in suitable coordinate sys- 
tem to compare with the data. Using the present theoretical model, we can calculate 



/ = 0. 



§4. Summary 
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observables in whole energy region of interest and can compare with the experimenal 
data directly. The final state interacton is partly included in the model since we have 
used the fit to the data of deuteron photodisintegration. 

Our model is found to reproduce gross features of all the existing data but 
certain descrepancies remain, which are (i) energy dependence of the total cross 
section for 3 He(7,pn)p sp and (ii) the step-like change of the total cross section of 
4 He(7,pn)ds P reaction around E 7 ~ 150 MeV. We found that the these discrepancies 
are difficult to reproduce by the present phenomenological model. Thus, we think 
that these discrepancies involve important information on the essential differences of 
the photodisintegration reactions of 3 He and 4 He from that of deuteron and should 
be treated carefully. 

For further studies, we need to develope a microscopic model based on the accu- 
rate wave functions of initial and final states and photon interaction with hadronic 
currents in order to understand these reactions deeper. We think it is very impor- 
tant to calculate the observables in the suitable form to compare with data by the 
microscopic model. 
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Figure Caption 

Figure 1: Double differntial cross sections of the 3 He(7,pn)p sp reactions for 
three differnt angles at four photon energies as indicated in figures. The 
solid and dashed lines indicate the calculated results withjijF=150MeV and 
p/=200MeV, respectively. Data are measured by TAGXlHP. 

Figure 2: Differntial cross sections of the 3 He(7,pn)p sp reactions at four 
photon energies as indicated in figuers. The solid lines indicate the calcu- 
lated results with pj=150MeV. Data are taken from ref.lilP. 

Figure 3: Total cross section of the 3 He(7,pn)p sp reaction is shown as a 
function of incident photon energy. The solid curveindicates the calculated 
result with py=150MeV. Data are taken from ref. tf . 

Figure 4: Total cross section of the 4 He(7,pn)d sp reaction is shown as a 
function of incident photon energy. The solid curve indicate^ the result 
with p/=200MeV. Data are those from refs. tSP (solid circles),^ (squares), 
(crosses) , © (open circles) , © (triangles) . 

Figure 5: Ratio of a( 4 He( / y,pn)d sp )/a(d( r y,p)n) as a function of E 7 . 
ai^He^ ,pri)d sp ) is calculated with p^=200MeV. Data are those from refs. 
(solid circles) , 0) (squares) , (crosses) , 0) (open circles) . 

Figure 6: Same as in Fig. 4. Each line indicates calculated result with 
p/=100(dashed curve), 200(solid curve), 300(dotted curve), and 400(dash- 
dotted curve)MeV, respectively. Calculated results are normalized to re- 
produce experimental peak height at A energy region. These normal- 
ization factors are 1.45 for p / =100MeV, 1.53 for p f 200MeV, 1.67 for 
p/=300MeV, and 1.88 for p /= 400MeV. 

Figure 7: Same as in Fig. 5. Each line corresponds to calculated results 
with different pf values and normalization factors as given in Fig. 6. 
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